close all;clear;clc;
% 中国各省、市、自治区经纬度(北纬、东经)
city=[31.52 117.17;% 安徽合肥
         39.55 116.24;% 北京
         29.59 106.54 ;% 重庆
         26.05 119.18;% 福建福州
         36.04 103.51;% 甘肃兰州
         23.08 113.14;% 广东广州
         22.48 108.19;% 广西南宁
         26.35 106.42;% 贵州贵阳
         20.02 110.20;% 海南海口
         38.02 114.30;% 河北石家庄
         34.46 113.40;% 河南郑州
         45.44 126.36;% 黑龙江哈尔滨
         30.35 114.17;% 湖北武汉
         28.12 112.59;% 湖南长沙
         43.54 125.19;% 吉林长春
         32.03 118.46;% 江苏南京
         28.40 115.55;% 江西南昌
         41.18 123.25;% 辽宁沈阳
         40.48 111.41;% 内蒙古呼和浩特
         38.27 106.16;% 宁夏银川
         36.38 101.48;% 青海西宁
         36.40 117.00;% 山东济南
         37.54 112.33;% 山西太原
         34.17 108.57;% 陕西西安
         31.14 121.29;% 上海
         30.40 104.04;% 四川成都
         39.02 117.12;% 天津
         29.39 91.08;% 西藏拉萨
         43.45 87.36;% 新疆乌鲁木齐
         25.04 102.42;% 云南昆明
         30.16 120.10;% 浙江杭州
         21.23 115.12;% 香港
         21.33 115.07;% 澳门
         25.03 121.30;% 台湾台北
         ];
save('./city_data.mat',"city",'-double');
a = 1:34;
a = randperm(length(a));
cityb=zeros(size(city));
for i =1:length(a)
    cityb(i,:) = city(a(i),:);
end
cityb =[cityb;cityb(1,:)];
figure('Position',[100 100 640 480])
hold on;
scatter(cityb(:,2),cityb(:,1),'r*');
plot(cityb(:,2),cityb(:,1),'b-','LineWidth',1.2);